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ABSTRACT 

We investigate the nature of gas accretion onto haloes and galaxies at z — 2 us- 
ing cosmological hydrodynamic simulations run with the moving mesh code AREPO. 
Implementing a Monte Carlo tracer particle scheme to determine the origin and ther- 
modynamic history of accreting gas, we make quantitative comparisons to an other- 
wise identical simulation run with the smoothed particle hydrodynamics (SPH) code 
GADGET-3. Contrasting these two numerical approaches, we find significant physi- 
cal differences in the thermodynamic history of accreted gas in massive haloes above 
~ 10 10 ' 5 Mq. In agreement with previous work, GADGET simulations show a cold frac- 
tion near unity for galaxies forming in massive haloes, implying that only a small 
percentage of accreted gas heats to an appreciable fraction of the virial temperature 
during accretion. The same galaxies in AREPO show a much lower cold fraction, for 
instance < 20% in haloes with Mhaio — 1O 11 M0. This results from a hot gas accre- 
tion rate which, at this same halo mass, is an order of magnitude larger than with 
GADGET, together with a cold accretion rate which is lower by a factor of two. These 
discrepancies increase for more massive systems, and we explain both trends in terms 
of numerical inaccuracies with the standard formulation of SPH. We note, however, 
that changes in the treatment of ISM physics - feedback, in particular - could modify 
the observed differences between codes as well as the relative importance of different 
accretion modes. We explore these differences by evaluating several ways of measuring 
a cold mode of accretion. As in previous work, the maximum past temperature of gas 
is compared to either a constant threshold value or some fraction of the virial tem- 
perature of each parent halo. We find that the relatively sharp transition from cold 
to hot mode dominated accretion at halo masses of ~ 10 11 M Q is a consequence of 
the constant temperature criterion, which can only separate virialised gas above some 
minimum halo mass. Examining the spatial distribution of accreting gas, we find that 
the filamentary geometry of accreting gas near the virial radius is a common feature in 
massive haloes above ~ 10 11-5 M©. Gas filaments in GADGET, however, tend to remain 
collimated and flow coherently to small radii, or artificially fragment and form a large 
number of purely numerical "blobs" . These same filamentary gas streams in AREPO 
show increased heating and disruption at 0.25-0.5 r v i r and contribute to the hot gas 
accretion rate in a manner distinct from classical cooling flows. 

Key words: cosmology: theory - galaxies: evolution, formation, haloes - methods: 
numerical 



1 INTRODUCTION 
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The manner in which forming galaxies acquire their baryonic 
matter through cosmic time in a ACDM cosmology is not yet 
fully understood. In the classic picture of accretion from the 
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intergalactic medium gas shock-heats to the virial tempera- 
ture of the halo and forms a hot pressure supported atmo- 
sphere in approximate equilibrium JRees & Ostriker 1977 



Silk 1977 White & Roes 1978). Gas can then cool onto a cen- 
trally forming galaxy, or not, based on the relevant cooling 
time-scale (White & Frcnk 19911. In this regime the accre- 
tion rate onto the galaxy depends on the cooling rate and 
not on the cosmological accretion rate onto the halo itself. 
However, this holds true only for sufficiently massive haloes. 
Below some threshold mass the dark matter halo cannot sup- 



port a stable hot, gaseous atmosphere (Birnboim & Dekel 



20031 and accreting gas does not necessarily shock heat to 
the virial temperature. This conclusion clearly motivates the 
differentiation of a "cold accretion mode" comprised of gas 
which does not shock heat to the virial temperature of its 
parent halo during accretion. Numerical simulations over the 
past decade have shown evidence for a cold accretion mode 
which dominates both at high redshift and for low mass ob- 



jects at late times i 


Katz et al.|2003 Keres et al.|2005 2009 


Ocvirk et al. 2008 


i. Haloes have been found to transition 



from cold to hot mode dominated at a mass in rough agree- 
ment with analytical estimates for the stable virial shock 
cutoff (e.g. |Birnboim fc Dekel|2003 l. 

Cold accretion in low mass objects below this thresh- 
old is precisely the regime where the virial shock is inac- 
tive. More massive haloes supporting a virial shock are also 
claimed to experience continuous, active cold mode accre- 
tion at high redshift [z ^ 2). This cold gas accretes along 
strongly filamentary geometries, termed "cold flows," which 
penetrate the virial radius and deliver gas deep within the 
halo ([ Keres et aL][2005l |Dekel et aL][2^09} |Dekel fc BlriT| 



boim 



2006) 



In the context of analytical predictions, the 
non-spherical accretion of cold gas via highly collimated ge- 
ometries was an unexpected result of these numerical simu- 
lations. Although smoothed particle hydrodynamics (SPH) 
simulations were the first to address the question of gas ac- 
cretion modes in a cosmological context ( Abadi et al.|200"3 



Keres et al. 2005 I they were soon followed by grid-based 



adaptive mesh refinement (AMR) efforts ( Ocvirk et al. 2008 
and others). 

Simulations of cosmological volumes using the SPH 
technique typically make a cold mode gas selection based on 
the past temperature history of each gas element, invoking 
the pseudo-Lagrangian nature of the scheme. The resulting 
maximum past temperature T max value is compared to some 
temperatu re threshold (|Keres et al 2005 , 2009 . van de Voort 



et al.|2011 Faucher-Giguere et al.|2011 I . This threshold was 
originally motivated by the location of a minimum in the 
strongly bimodal global distribution of T max found by|Kere§| 



et al.| 1 2005 1, using the TreeSPH ( |Hernquist fc Katz||l989| ) 
code. That work compared T max to the virial temperature 
of haloes as well as to a constant temperature threshold, and 
chose T c ~ 2.5 x 10 5 K as sufficient to identify virial heating 
in haloes, due to the high halo masses considered. Essen- 
tially all SPH simulations to date have used this approach, 
although with less care as to its applicability, as we discuss 
later. One notable exception is |Brooks et al.| (|2009[| w hich 
used the SPH code GASOLINE ( |Wadsley et al.||2004j > and 
identified shocked gas particles based on an entropy jump 
criterion. This is a novel approach which should be explored 
in a full cosmological volume. The authors concluded that 
their method led to a selection of shocked/unshocked gas 



which agreed well with the selection of hot/cold gas based 
on the use of a constant temperature threshold. 

More recently, finite volume techniques including cos- 
mological re-simulations using adaptive mesh refinement 
(AMR) have been used to study the same problem. Grid- 
based simulations offer a better treatment of shocks, fluid 
instabilities, and phase boundaries. However for the analy- 
sis of cosmological gas accretion they are limited by their 
Eulerian nature, as the trajectory of any gas parcel cannot 
be simply identified due to mass exchange between cells. 
This difficulty can be overcome using tracer particles, and 
we argue these are a necessary feature of modern cosmo- 



logical simulations. One AMR work to date (Dubois et al 



2012 1 used tracer particles in order to study the angular mo- 
mentum of accreting gas. A common alternate approach is 
to use the instantaneous gas properties such as mass flux 
through a surface (e.g. |Ocvirk et al.|2008"l|Dekel et al.|2009 1 
to measure rates, or the low metallicity of gas as an indica- 



tor of cosmological origin ( Agertz et al. 2009 1. Unfortunately 



these techniques cannot constrain the past thermal history 
of any gas element, nor identify for example material tidally 
stripped from a galaxy or subhalo during a previous strong 
interaction. Given these caveats, however, AMR simulations 
of single haloes have also found strong support for the im- 
portance of a filamentary, cold accretion mode in massive 
systems at high redshift. 

The contribution of "clumpy" accretion due to infalling 
satellites and other bound substructures is also difficult to 
distinguish in Eulerian simulations. Some SPH studies use a 
stringent definition based on lack of membership to any iden- 
tified bound substructure at any time prior to accretion (e.g. 
|Brooks et al.|2009] ). Others only require non-membership at 



some particular previous time (Keres et al. 2005 van de 



Voort et al. 2011). In either case, this gas should not be 
considered as contributing to cold accretion from the inter- 
galactic medium (IGM). This is appropriate since we wish 
to isolate cosmological gas accretion from the IGM and not 
hierarchical buildup via mergers. The technical distinction is 
not particularly important as smooth accretion (defined as 
all gas not identified as clumpy) is generally found to be the 
dominant contribution to the total accretion rate ( |van de| 



Voort et al. 2011 1. However, Shen et al. (20121 for example 



includes a substantial ~35% contribution of infalling dwarf 
satellites to the cold accretion mass budget. This effectively 
combines gas from multiple sources under a single label and 
ignores the diverse physical processes driving distinct accre- 
tion mechanisms. This is especially important in the results 
we report in what follows using AREPO, for which the ac- 
cretion rate of cold gas directly from the IGM is small. 

Only recently have simulations begun to explore the 
effects of feedback and galactic scale outflows on the charac- 
ter of gas inflow ( Oppenhcime r et al.|2010| |Faucher-Giguere| 
et al.|2011 Ivan de Voort et al. 120111). Differences in morphol- 



ogy between different modes of accretion in massive haloes 
and the added complexity from the physics of star formation 
implies that accretion rates may be quite sensitive to feed- 
back processes. Simulations using SPH have concluded that 
a filamentary cold mode is largely unaffected by the presence 
of strong outflows arising from blastwave supernova feed- 
back ( |Shen et aL"1|2012[ ). |van de Voort et al.| ( |2011[ ) studied 
the impact of AGN feedback on inflow and concluded that 
it preferentially prevented hot mode gas (under the stan- 
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dard definition) from cooling from the halo onto the galaxy, 



while Dubois et al. (2012) found that AGN feedback mor- 



phologically disturbed cold filaments. Current treatments of 
feedback in cosmological simulations, however, are largely 
phenomenological and are not self-consistent. Furthermore, 
the same numerical issues explored in Vogelsb erger et ah] 
(2012]), |Sijacki et aT] pfJ12|, |Kere§ et al.| (|2012[> and dis- 



cussed herein that compromise the accuracy of SPH studies 
of gas inflow will also affect the interaction of outflowing 
ejecta and wind material with both halo and filamentary 
gas. Moreover, the accretion rates onto galaxies are signif- 
icantly higher in simulations with AREPO than with GAD- 
GET, implying that feedback effects need to be even more 
efficient than in previous simulations with SPH. 

Strong feedback may also indirectly affect the charac- 
terisation of inflows. In particular, galactic fountain mate- 
rial and other recycled gas may be a significant additional 



mode of accretion onto galaxies. Oppenheimer et al. ( 2010 1 



included a phenomenological galactic wind model in cos- 
mological SPH simulations and concluded that recycled gas 
accretion is a prominent "third" mode which in fact is the 
dominant accretion mechanism at z ^ 1, with minimal ef- 
fect on this high redshift accretion with choice of outflow 
velocity. The ejection of hot gas from galaxies might also 
physically be expected to add mass and thermal support to 
a halo atmosphere. The result would be increased heating of 
cosmologically inflowing gas, an underestimated effect in all 
simulations lacking galactic scale outflows, including those 
presented here. 

This paper presents results from a comparison of cos- 
mological simulations run with the SPH code GADGET 



(| Springel 2005 1 and the moving mesh code AREPO ( Springel 



2010 1. We compare individual objects as well as statistics de- 



rived from cosmological volumes and show that the impor- 
tance and character of cold gas accretion depends sensitively 
on both analysis methodology and numerical technique. In 
Section[2]we describe the numerical details of our approach, 
two complementary tracer particle schemes, and our analy- 
sis methodology. Section[3]compares the thermodynamics of 
gas accretion between the two codes, discusses the various 
definitions of the cold mode and the sensitivity of our con- 
clusions to those definitions, and investigates the geometry 
and filamentary nature of inflowing gas. Finally, Section [4] 
wraps up our discussion and summarizes the main conclu- 
sions. 



2 METHODS 

We use and extend the simulations and numerical methods 



presented extensively in Vogelsberger et al. (20121 and refer 



the reader to that paper for additional details. Brief descrip- 
tions of the key hydrodynamic differences between the two 
codes are presented here. 

GADGET (last described in |Springel| |2005[ ) is a 
smoothed particle hydrodynamics (SPH) implementation 
where the Euler equations are solved to model fluid flow by 
discretising mass using a set of particles. Continuous fluid 
quantities are defined by an interpolation kernel over a set of 
nearest neighbors. This method is pseudo-Lagrangian in na- 
ture (see discussion) and is particularly well-suited to handle 
the large dynamic range present in cosmological simulations 



where it is desirable to resolve both large scale structure 
and the internal dynamics of forming galaxies. We use a 
"standard" SPH formulation in GADGET to facilitate com- 
parison to earlier work ( jMona ghan 1992J) - the cubic spline 
kernel with 32 neighbors and the density based formulation 



(Springel & Hernquist 2002) modified via an artificial vis- 
cosity osph = 1-0 with the Balsara switch ( Balsara|1995 I. 

AREPO (described in Springel 2010) is a finite volume 
scheme where the control volumes are defined by a Voronoi 
tessellation of space. Euler's equations are solved using Go- 
dunov's method with the MUSCL-Hancock scheme to com- 
pute numerical fluxes and obtain second order accuracy. The 
tessellation is obtained by a set of mesh generating points 
which are allowed to move arbitrarily, though in our simu- 
lations follow the flow in a quasi-Lagrangian fashion. As a 
result AREPO retains the principal strengths of SPH includ- 
ing its adaptivity to a large dynamic range in spatial scales, 
Galilean invariance of the truncation error, and an accurate 
and efficient gravity solver (e.g. O'Shea et al.|[2005 k It also 
gains the strengths of finite volume codes, including the im- 
proved treatment of fluid instabilities, weak shocks (which 
can be missed in SPH; (e.g. |Keshet et al. 2003)), phase in- 
terfaces, and shearing flows. We use the maximum cell face 
angle regularisation scheme and allow dynamic refinement 
and de-refinement to maintain approximately constant mass 
cells ( Vogelsberger et al.|2012 1 . 



2.1 Simulation Set 

All our simulations employ a WMAP-7 cosmology (Q.a,o = 
0.73, fi m , = 0.27, tt bfi = 0.045, a 8 = 0.8 and h = 0.7) with 
a 20h~ 1 Mpc sidelength volume at several different numeri- 
cal resolutions, initially with equal numbers of dark matter 
particles and gas elements. The particle counts, initial gas 
element masses, dark matter particle masses, and Plummer 
equivalent gravitational softening lengths are given in Table 
|TJ. Unless otherwise noted, all results in this paper com- 
pare results between the 512 3 resolution runs at a redshift 
of z = 2. 

GADGET and AREPO simulations include identical 
physics. We account for optically thin radiative cooling as- 
suming a primordial H/He ratio which sets an effective tem- 
perature floor for gas in these simulations just below ~ 10 4 



K (Katz et al. 1996 1. We do not include metal line cooling. A 



redshift-dependent ionizing UV background field (Faucher- 
|Giguere et al.|20 09) is included as a spatially uniform heat- 
ing source. Star formation and the associated ISM pressuri- 
sation from unresolved feedback events are included follow- 
ing [SprjngeT&Tlernquist (2003]). Gas elements are stochas- 
tically converted into star particles of a constant mass when 
the local gas density exceeds a threshold value of riu ~ 0.13 
cm" 3 . 

We do not include any explicit forms of strong stellar 
feedback that would drive galactic-scale winds. There is no 
treatment of black holes, radiative transfer and its effects, 
or magnetic fields. The key point is that we use a "standard 
physics" implementation which has been well studied and 
probes the interaction of gravity, hydrodynamics, and star 
formation in the cosmological context. The gravity imple- 
mentation as well as the chosen subgrid physics are identical 
between our SPH and moving mesh simulations, allowing us 
to attribute differences in the solutions to differences in the 
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Table 1. Details on the three different simulation resolutions and relevant numerical parameters. 



Gas Elements 


DM Particles 


Vel Tracers 


MC Tracers 


m targct/SPH [ h lM @] 


m DM [h 1 M ] 


e [h- 1 kpc] 


128 3 


128 3 


1 x 128 3 


10 x 128 3 


4.8 x 10 7 


2.4 x 10 s 


4.0 


256 3 


256 3 


1 x 256 3 


10 x 256 3 


6.0 x 10 6 


3.0 x 10 7 


2.0 


512 3 


512 3 


1 x 512 3 


10 x 512 3 


7.4 x 10 5 


3.7 x 10 6 


1.0 



method used to solve the equations of hydrodynamics. They 



are also similar to Keres et al. (20091, in terms of the code 



and physics employed as well as the simulation resolution, 
enabling us to make a direct comparison to that work. 



2.2 Tracer Particles 

In order to trace the evolution of gas properties over time in 
our moving mesh calculations we use a new "Monte Carlo" 
tracer particle technique (see Vogelsberger et al. and Genel 
et al., in prep). This probabilistic method associates tracers 
with parent gas cells and exchanges them based explicitly 
on the mass fluxes through each face. Specifically, all cells 
start with an equal number of tracers, after which each child 
tracer transfers from its original cell (i) to a neighboring cell 
(j) with a probability equal to AMy/M;, the ratio of the 
mass flux through face (ij) to the current mass of the origi- 
nating cell. By construction, the tracer density is guaranteed 
to follow the underlying fluid density, at the cost of Poisson 
noise due to the probabilistic nature of the scheme. Further- 
more, at each computational timestep each tracer records 
several fluid quantities of its parent gas cell - for this work 
we use only the maximum previous temperature and time of 
that Tmax event. We initialize each gas cell with 10 Monte 
Carlo tracers. Due to the target mass refinement scheme 
the number of child tracers per parent cell remains roughly 
constant throughout the simulation, with a Gaussian distri- 
bution extending a factor of two above and below the initial 
number. In particular, in galaxies at z — 2 the mean number 
of child tracers is ~10 with a standard deviation of ~3. 

We also evaluated a "velocity field" tracer particle 
scheme (see Vogelsberger et al. and Genel et al., in prep). In 



this more typical approach (e.g. Vazza et al. 2010 Seiten- 
zahl et al~|2 010 Du bois et al.|2012 |, tracers represent mass- 
less, passive particles which are purely advected by the local 
velocity field of the fluid. Our testing shows that this ap- 
proach exhibits a systematic bias in its Lagrangian ability 
to follow mass fluxes. We suggest that the technique does 
not accurately recover the flow of mass in astrophysical sit- 
uations involving convergent flows, including cosmological 
cases where the problem manifests as a "pile-up" of trac- 
ers in the centers of dark matter haloes (see Genel et al., 
in prep). The Monte Carlo approach is unaffected by this 
bias seen with the velocity field tracer particles. Further- 
more, the flux-minimising moving mesh scheme allows us to 
use a reasonably small number of such tracers while achiev- 
ing an acceptable level of noise (Genel et al., in prep). All 
results presented herein which involve tracking gas proper- 
ties through time in the AREPO runs use our Monte Carlo 
tracer method. We however checked that our conclusions are 
qualitatively robust against choice of tracer scheme. 



2.3 Identifying Haloes 

We identify dark matter haloes using the SUBFIND algo- 
rithm (Springel et al. :2001) which begins with a friends-of- 
friends procedure (linking length b — 0.2) applied to dark 
matter particles. Gas and stars are associated to their near- 
est dark matter particle. Tracers are naturally associated 
with their nearest (parent) gas cells. Next, an iterative un- 
binding procedure which accounts for thermal energy identi- 
fies substructures within each FOF group which are gravita- 
tionally bound groups, with a minimum of 20 particles each. 
We refer to the largest such subgroup as the halo itself. Thus 
the halo selection specifically excludes gas in satellites and 
other orbiting substructures as well as unbound material, 
at the time of selection. However, it is important to note 
that material associated but not bound to substructures, 
including tidal features such as tails, would not be explicitly 
excluded by this procedure. We discuss below the separate 
issue of classifying gas as either smooth or clumpy when 
following its thermal history. 

The center of each halo is taken as the position of the 
particle with the minimum gravitational potential. We take 
as the virial radius r v i r the numerically computed r2oo,cri<F] 
and for the halo mass the mass of the largest bound sub- 
group. To determine accretion times, we need to estimate 
the time evolution of the virial radius of each halo as well 
as its position. However, SUBFIND is a purely 3D algorithm 
with no explicit time linking. We construct a parent tree (a 
simplified merger tree) by matching haloes between succes- 
sive snapshots based on the largest fraction of common dark 
matter particles. We require 60% particle agreement and 
reasonable mass and center offsets between snapshots. For 
consistency we restrict our analysis on a halo by halo basis 
to the time period for which the parent tree is well denned, 
which is true for > 90% of haloes with M ha i > 10 10 M Q over 
the chosen time window, as discussed below. 



2.4 Maximum Past Temperature 

We broadly consider three different means of measuring cold 
mode accretion. Two depend on the properties of associated 
haloes. We introduce three important temperature values 
used throughout this work: 

• Tvir.cur is the virial temperature for the parent halo of 
each tracer or gas particle at the "current" time (i.e. the 
time the gas selection is made, z = 2). 

• 3vir,acc is the virial temperature evaluated at the "ac- 
cretion" time on a particle by particle or tracer by tracer 



1 The radius within which the enclosed overdensity is 200 times 
the critical density p cr it. 
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Figure 1. The mass- weighted histogram of past maximum temperature T max for gas accreted onto central galaxies (top) and halo 
atmospheres (bottom) by z = 2, as a function of the parent halo mass at z = 2. AREPO and GADGET are shown separately in the left 
and right panels, respectively. In this figure we do not normalize these temperatures by any fraction of T v i r cur or T v i r acc . A constant 
temperature of T c = 10 5 ' 5 K is shown (dashed line) together with the virial temperature at redshift two (solid line). The two codes 
deviate strongly for galaxies in haloes with masses above ~ 10 10,5 Mq, where the strong contribution of hot gas in AREPO, which scales 



approximately with T v j r , is mostly absent in GADGET. The cold 
largely absent in the AREPO simulation. 



10 K gas which dominates the accretion in GADGET galaxies is 



basis. We define the accretion time here as the most recent 
virial radius crossing time. 

• Tc refers to a constant temperature threshold at some 
value, applied uniformly across all halo masses. 

To estimate the virial temperature we use the definition 



of Barkana & Loeb (2001 1 which at redshift two gives 



T ■ — 

-L vir — 



2k B 



4 x 10 b K 



M, 



10 11 h^M 



2/3 



(1) 



where V c is the circular velocity at the virial radius and 
H ~ 0.6 for fully ionized, primordial gas. We note that the 
redshift scaling as expressed here is only approximate, as 
it ignores small changes in cosmological parameters at the 
redshifts of interest. 



We calculate a cold fraction for each halo or galaxy as 
the ratio of cold to hot accreted gas at z = 2 using each 
of the three methods. To calculate the maximum past tem- 
perature of a gas selection in GADGET runs we examine 
the temperature of each SPH particle at each snapshot back 
to the start of the simulation. We use 314 snapshots loga- 
rithmically spaced in redshift from z = 30 to z = 0, cor- 
responding to a mean spacing at z — 2 of ~ 25 Myr. To 
calculate the maximum past temperature of a gas selection 
using the Monte Carlo tracers in AREPO runs we first find 
all the tracer children of those gas cells at the target redshift. 
These tracers are then treated as a uniform population and 
the maximum past temperature of each, recorded for each 
active timestep of their parent gas cell, is obtained since the 
start of the simulation. That is, the mass weight of the cell 
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is divided among its children tracers and each may poten- 
tially have a distinct thermodynamic history. In both cases 
we neglect temperature maxima while gas is on the effective 
equation of state (star forming). 



2.5 Measuring Accretion 

We separate gas in the "central galaxy" from the "halo at- 
mosphere" by a cut in the density, temperature plane as 

1 



log(T gas [K] 



■ log (p gas [10 iu M Q fe'kpc" 3 ]) < 6.0 



(2) 



following |Torrey et al.| ((201 1 ) . This selection differentiates 
between hot gas in the halo and the cold, dense gas which is 
rotationally supported at the center of the halo. We add a 
radial cut at 0.15r v i r to this separation criterion to prevent 
including any cold gas at large radii as part of the galaxy. 
Star particles are also included throughout our analysis and 
are considered as part of the galaxy if they satisfy the radial 
cut. Our results are quantitatively unaffected for reasonable 
choice of these three galaxy selection parameters. 

In order to define accretion and measure accretion rates 
we introduce a "time window" over which accretion events 
are counted. A gas element is considered to have accreted 
onto a halo if it belonged to that halo at z — 2 and crossed 
the virial radius during the chosen time window. A gas ele- 
ment is considered to have accreted onto a galaxy if it be- 
longed to that galaxy at z = 2 and either crossed the phase 
space cut in (p,T) or the radial cut at 0.15r v i r during the 
chosen time window. The accretion time is taken to be the 
more recent (lower redshift) of these two events. The same 
definition is used for star particles in galaxies, where we 
search for the accretion event of the progenitor of each star 
particle only in the gas phase. That is, stars which remain 
stars for the duration of the chosen time window are not 
considered as accreted. Rates are calculated by normalising 
the total accreted mass (counting particles or tracers) by the 
time period. 

All accreted material is separated into one of three dis- 
joint "modes" of accretion based on the following criteria, 
evaluated at the virial radius crossing time for each gas ele- 
ment: 

• Smooth - not a member of any halo or subhalo, other 
than the main progenitor branch of its parent halo, and like- 
wise back to z — 6. 

• Clumpy - gravitationally bound to a halo or subhalo, 
other than the main progenitor branch of its parent halo. 

• Stripped - not a member of any halo or subhalo 
(smooth), but gravitationally bound to some other halo or 
subhalo at any previous time back to z = 6. 

These definitions allow us to exclusively consider accre- 
tion from the IGM by restricting our analysis to the smooth 
mode only. This is required in order to remove the "merger 
contribution" , and is similar in spirit to the approach taken 
by previous studies of cosmological gas accretion for which it 
is possible to make this distinction (SPH and not grid codes). 
However, we note that our definition of clumpy is somewhat 



more restrictive than in previous works (e.g. Keres et al. 
2005) which would remove only the central galaxy compo- 



nent of each subhalo. By taking all gas bound to substruc- 
tures as a merger contribution we impose a stronger restric- 



tion, and therefore remove a larger fraction of material in 
our analysis of smooth accretion. Finally, we measure the 
importance of tidally stripped material with the third mode, 
which is not direct accretion from the IGM but rather gas 
which has been previously gravitationally bound to some 
halo or galaxy. Unless otherwise noted, all results in this 
paper include only smoothly accreted material over an ac- 
cretion time window of 1 Gyr. The only exceptions are Fig- 
ure [6] which considers the total accreted baryonic mass over 
an "integrated" time window extending back to z — 6, and 
Figure [11] which compares the smooth, clumpy, and stripped 
accretion modes. 



3 RESULTS 

In this section we investigate differences between GADGET 
and AREPO in the thermodynamic history of gas accreted 
onto haloes and the galaxies forming at their centers. Our 
primary proxy for this potentially complicated history is 
Tmax, the maximum past temperature of each gas element. 
In Figure [l] we show the mass- weighted histogram of T max 
for smoothly accreted g a function of parent halo mass. 
We also indicate a nominal constant temperature threshold 
value along with the virial temperature at redshift two. 

Firstly, we see that GADGET and AREPO agree well on 
the relative distribution of past maximum temperatures con- 
tributing to halo atmospheres across the entire mass range 
(lower panels). For accretion onto galaxies (upper panels) at 
all halo masses the AREPO temperature distribution scales 
approximately with T v i r , while in GADGET this behavior 
holds only for low mass haloes below ~ 1O 1O M0. In fact, we 
find that gas heating proceeds in relation to the virial tem- 
perature even for low mass systems. At ~ 1O 9 M0 the T max 
distribution may be beginning to scale somewhat shallower 
than Tvir, flattening towards a uniform, cold temperature 
set by the IGM of a few times 10 4 K. However, our simula- 
tion resolution is not sufficient at these low mass systems to 
convincingly demonstrate this. There is also an indication of 
a roughly constant temperature component at T max ~ 10 5 K 
in the GADGET simulation, which we discuss further below. 
The scaling with T V - 1T , together with the lack of any distin- 
guishable segregating feature about a constant temperature 
threshold on the order of a few times 10 5 K, motivates the 
comparison of T max to T v i r on a halo by halo basis. 

In Figure [2] we show the distribution of T max for 
smoothly accreted gas, onto both galaxies and haloes, di- 
vided into six bins of halo mass. In each case we now normal- 
ize temperatures by the virial temperature of the dark mat- 
ter halo of each gas element, at the time of accretion. For ac- 
cretion onto central galaxies for haloes with M 1O 1O ' 5 M0 
the two codes show excellent agreement. For such low mass 
systems, however, this agreement does not necessarily indi- 
cate a similar level of shock heating. The virial temperatures 
of these haloes are sufficiently low that IGM pre-heating be- 
gins to set an effective temperature floor for the least mas- 
sive systems. This heating outside the halo occurs through 
a combination of radiative energy input from the UV back- 
ground, as well as shocks during large scale collapse, both 
of which we discuss further in the context of cold fractions. 
It is only for Mh a i ^ 1O 1O ' 5 M0 where gas accreted onto 
centrally forming galaxies begins to show different thermal 
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Figure 2. Relative distributions of past maximum gas temperature normalized by the virial temperature of the parent halo at the time 
of accretion, without any cold or hot mode differentiation. Each of the four curves in each panel is separately normalized, so that only 
their positions and shapes, but not their amplitudes, should be compared. At low halo masses (top three panels) and for gas accreted onto 



haloes (dotted lines) the two codes agree well. For more massive haloes above 



10 1 



5 Mq significant differences arise in the thermal 



history of gas accreted onto galaxies, where gas in GADGET has predominantly T max < T v ; r acc while gas in AREPO has predominantly 
?inax ^ T v ; r ace- 



histories between GADGET and AREPO. At these masses, 
the majority of gas in the SPH calculation is never heated to 
temperatures comparable to the halo virial temperature at 
accretion. In contrast, a dominant fraction of the same ac- 
creted gas in AREPO reaches approximately the virial tem- 
perature, or above. This striking difference in the temper- 
ature history of gas accreted onto galaxies in haloes with 
Afhaio 10 10 ' 5 Mq is the first main result of this paper. 

In this same mass regime both codes show evidence 
for potentially bimodal temperature distributions. This is 
most evident for GADGET galaxies between 10.5 ^ log(M) 
^ 11.0 where T max /T v i r ,acc could be well described by the 
sum of two symmetric distributions with different peak tem- 
peratures. As normalized, the colder peak moves to lower 
temperatures with increasing halo masses approximately as 
T oc M 2,/3 and is consistent with a constant physical temper- 
ature of T ~ 10 5 K, as seen in Figure [l] We note that com- 
mon choices for a constant temperature threshold adopted 
by previous studies are above this value. The warmer peak 
shows no scaling with halo mass and is consistent with a 
heating process scaling with virial temperature. These two 
contributions represent, respectively, accretion of roughly 
constant temperature material from the IGM which has min- 
imal interaction with the halo potential, and material which 
virialises to a mass dependent temperature. Results from 
both GADGET and AREPO are consistent with this picture, 



though they strongly disagree on the relative contributions 
of these two accretion channels. 



3.1 Accretion Rates 

To understand the nature of the differences, we measure the 
accretion rate onto both central galaxies and haloes them- 
selves as a function of halo mass. Figures [3] and [4] sepa- 
rate cold from hot mode accretion based on comparison of 
Tmax to either a constant temperature threshold T c or the 
halo virial temperature T v ir, respectively. The lower panels 
of both figures again indicate that the total accretion rates 
and thermal history of gas which accretes onto the halo it- 
self is similar between the two codes. Our accretion rates 
for the low mass end are somewhat lower than in IKeresI 



et al. ( 2009 1 for instance, which may be due to our more 



stringent definition of smooth accretion as well as our bet- 
ter resolution of the hierarchical nature of structure growth. 
When using a Tmax/Tvir.acc comparison this scaling is well 
described by a power-law in the mass range where systems 
are both well resolved and sufficiently numerous. We now 
discuss how the accretion rates compare between the two 
codes based on Figure [4] and T v i r , a cc, returning later to this 
particular choice. 

The accretion rate onto galaxies in GADGET is domi- 
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Figure 3. Gas accretion rate onto galaxies (top) and haloes (bottom). In all four panels we define the cold mode as requiring T max be 
less than a constant value T c , as listed, demonstrating the sensitivity of the derived rates to the method of measuring the cold versus 
hot mode. 
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Figure 4. As in the previous figure, exploring the sensitivity of the gas accretion rate onto galaxies (top) and haloes (bottom) to 
definition of the cold mode. Here we take cold mode gas as having T ma x/T v i riacc < 1.0,0.8,0.4, as listed. 
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Figure 5. Sensitivity of the derived cold fraction to the manner in which it is measured. The top two panels consider gas accreted onto 
central galaxies whereas the bottom two panels consider gas accreted onto halo atmospheres. The left two panels consider a gas element 
to have been accreted cold if the fraction T max /T v ; r acc is less than {1.0,0.8,0.4}. The right two panels define cold mode gas as having 
log (T max [K]) less than a constant value of {5.7, 5.5, 5.3}. We find that the most massive systems which have / co id — 1 in GADGET have 
a cold fraction approaching zero in AREPO, while the relatively sharp transition from cold to hot mode dominated, seen clearly in the 
constant temperature moving mesh case, is not evident when using the virial temperature comparison. 



nated by gas which has always been relatively cold, across 
the entire mass range of 10 10 < Mhaio < 10 12 . In AREPO the 
same galaxies are predominantly fed by gas which has ex- 
perienced significant heating. This contrast is robust under 
any reasonable comparison of maximum past temperature to 
either T c or T v j r . The accretion rate in the hot mode is dif- 
ferent by roughly an order of magnitude at Mh a i ~ 1O 11 M0 
between the two codes, having essentially zero contribution 
to haloes above this mass in the SPH simulation. The pri- 
mary reason for this difference is the efficient cooling of hot 
halo gas in AREPO which is artificially suppressed in GAD- 



GET by spurious viscous heating. Bauer & Springel (20121 



studied the dissipation of energy from subsonic turbulence 
in SPH and found that spurious heating arises both from 
the premature dissipation of turbulent energy on large spa- 
tial scales and the viscous damping of SPH-noise on small 
scales. In the context of our current cosmological simula- 
tions, this prevents gas at the cooling radius from forming a 
proper Kolmogorov-like spectrum. The viscous heating off- 
sets part of the gas cooling, keeping the gas locked in the 
diffuse, hot phase. In AREPO this same turbulent energy 
is more realistically dissipated by cascading to smaller spa- 
tial scales and higher densities, allowing significantly more 
efficient cooling of hot gas in the atmospheres of massive 



haloes (Keres et al. 20121. This cooling channel is a domi- 



nant accretion mode of massive haloes in the moving mesh 
calculation. 

The accretion rate of cold gas onto galaxies also differs 



between the two calculations. At M ha io 10 11 M© GADGET 
galaxies have a factor of two larger cold accretion rates as 
corresponding AREPO galaxies, and this discrepancy grows 
with increasing halo mass to an order of magnitude effect for 
Mh a io — 10 12 . We identify two primary numerical issues that 
cause the observed artificial increase of the cold gas accre- 
tion rate in the SPH calculation. Both operate by increasing 
the efficiency with which gas is heated in various regimes, 
resulting in a redistribution of gas to higher temperatures 
in the AREPO haloes. 



First, Creasey et al. (20111 demonstrate the issue of 



numerical overcooling due to shock broadening which arises 
due to the artificial viscosity treatment in SPH. At com- 
parable resolution, the shock-capturing Godunov scheme in 
AREPO better resolves the resulting post-shock gas proper- 
ties (including temperature), critical since smoothly accret- 
ing gas undergoes multiple small shocks as it virialises with 
hot halo gas. Secondly, spurious pressure forces in regions 
of steep density gradients ( Agertz et al.|2007 l, particularly 
contact discontinuities, lead to the numerical fragmentation 
of gas structures into artificial "blobs" due to surface ten- 
sion (see Torrey et al.|2011 i. A discussion of the formation 
mechanism of these blobs and their prevalence in different 
formulations of SPH can be found in |Hobbs et al.| ( |2012[ > . 
Their presence is a significant contribution to the accretion 
rate of cold gas in GADGET. We have checked explicitly 
that gas in these blobs have sufficiently low T max to be in- 
cluded in the cold accretion rate regardless of definition. 
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For Mhaio ^10 Mq these radially penetrating blobs con- 
tribute to the cold, smooth accretion rate in GADGET, di- 
rectly transporting low temperature, low angular momen- 
tum material into the centers of haloes (e.g. Keres & Hern- 
quist||2009 1. No comparable gas features are present in the 
moving mesh calculation, and this cold gas contribution to 
galactic accretion is entirely absent. 

3.2 Cold Fractions and the Transition Mass 

The ratio of total accreted mass with T max below some tem- 
perature threshold to the total mass above - the cold frac- 
tion, as a function of halo mass - is a widely reported quan- 
tity in studies of cosmological gas accretion, dating back to 
Keres et al. ( 2005 1. Figure[5]explores how robust the derived 
cold fractions are to different types of parameter choices in 
the measurement of a cold mode. Using a constant temper- 
ature comparison (upper right panel), we note that the cold 
fraction for galaxies derived from our GADGET simulations 



10 10 to 10 12 ; 



Mq. Galaxies in low mass haloes 



are in agreement with Keres et al. (20091. As expected from 
our previous discussion of accretion rates, this cold frac- 
tion is near unity for all halo masses, while in AREPO the 
cold fraction drops to zero for galaxies hosted in sufficiently 
massive haloes above ~ W 11 Mq. We have seen already that 
lower cold fractions in AREPO arise from both a lower total 
gas mass accreted through a cold channel, as well as a larger 
total gas mass accreted through a hot channel. 

We define the "transition mass" as the halo mass at 
which the cold fraction equals 1/2. Using the constant tem- 
perature cut T c ~ 5.5 ( |Keres et al.|2005 l leads naturally to 
a transition from / co id = 1 to / co id = at a halo mass of 
~ 10 110 M o in our highest resolution simulation. It is tempt- 
ing to compare this to model predictions for the critical halo 
supporting a virial shock. |Birnboim fc D ekcl ( 2003]) showed 
that this transition mass, defined as the point where the cold 
and hot modes have equal contribution, scales due to the 
contribution of metals to the cooling curve from ~ 1O 1O ' 9 M0 
at Z = to ~ 10 n - 6 M s at Z = 0.3 (redshift two). How- 
ever, we observe that the appearance and location of this 
transition is by construction due to the choice of T c . Given 
the scenario where all haloes heat some large fraction of in- 
falling gas to their virial temperature, this use of a constant 
temperature threshold would also reproduce a transition in 
this mass range due solely to the changing virial tempera- 
tures. Consider the case of a halo with a virial temperature 
just below the chosen constant temperature threshold, cor- 
responding to a mass of ~10 n MQ taking the usual choice 
for T c . Even if this halo shock heats all accreting gas to 
its virial temperature it would nevertheless be attributed a 
cold fraction of 100% under this method of measuring the 
cold accretion, which is clearly not a physically meaningful 
statement. We find that the subsequently derived transi- 
tion mass increases monotonically with the chosen T c if a 
constant temperature criterion is used to separate the two 
modes. From T c = 5.3 to T c = 5.7 (a factor of ~ 2.5) the 
transition mass increases by roughly half a decade, which is 
consistent with the y 3 / 2 scaling expected solely due to the 
increasing virial temperature. 

If we instead use a virial temperature comparison (up- 
per left panel), the derived cold fractions show markedly 
different behavior. In this case GADGET galaxies exhibit a 
monotonic increase in their cold fraction from ~ 0.4 to 1.0 



from M ha i 

below Afhaio ^ 10 1U M Q are no longer cold mode dominated, 
while those in high mass systems are. In contrast, AREPO 
galaxies still show the broad trend of decreasing cold fraction 
with increasing halo mass, though it tends towards contin- 
ual evolution as opposed to sharp transition. Using the virial 
temperature as the dividing criterion on a halo by halo ba- 
sis makes any transition significantly broader in halo mass. 
The second key result of this paper is then, given this more 
physically motivated means of measuring cold mode accre- 
tion in lower mass haloes, there is no longer a clear feature 
in the cold fraction as a function of halo mass. 

Although we argue for a comparison against T v h- and 
not T c in this work, it is important to note that the two 
methods work equally well within a specific halo mass range. 
In particular, for haloes with T v i r 3> T c such that virialisa- 
tion related heating is well separated by this cut. At z = 2 
this must require as a minimum Mh a i ^ 10 Mq which 
is indeed the mass range considered in early studies (e.g. 
Keres et al.||2005 l. Unfortunately this minimum halo mass 
falls well within the cold-hot transition range as reported at 
this redshift as well as at redshift zero, where the restriction 
becomes Mh a io ^ 10 12 Mq - the midpoint or even tail of the 
transition (van de Voort et al.||2011 Faucher-Giguere et al 



20111. Caution must also be taken using any T v i r compari- 



son as the more restrictive selections of the cold gas selection 
(0.4, 0.8) for low mass haloes fall below the effective temper- 
ature floor of our simulations. We run a separate experiment 
at the 256 3 resolution which excludes the heating contribu- 
tion from the UV background. The resulting temperature 
distribution of IGM gas extends to significantly lower tem- 
peratures, and the turnover in cold fraction around halo 
masses of ~10 10 ' 5 Mq essentially disappears. We then at- 
tribute this feature in the cold fraction to the high IGM gas 
temperature prior to accretion. This also indicates the mass 
scale below which a temperature criterion of any kind will 
be unable to separate physically distinct accretion modes. 

Finally, we note that switching from the "current" (not 
shown) to the "accretion" virial temperature criterion leads 
to a ~ 10% change, lowering the cold fraction across all 
halo masses. This change is unlikely to have great impact 
on the general interpretation of the cold mode importance. 
However, the virial temperature of a particular halo can 
change significantly between the time material is accreted 
and z = 2. To physically motivate our analysis requires that 
we either consider the virial temperature when it was rel- 
evant to the question of shock heating - at the time of ac- 
cretion - or restrict the analysis window to a short enough 
time-scale over which the virial temperature does not appre- 
ciably change. 



3.3 Integrated Baryonic Budget 

The strongly discrepant history of gas which eventually 
builds up galaxies by z = 2 has important ramifications for 
the morphology and structure of these forming galaxies. In 
Figure [6] we show the total gas mass smoothly accreted into 
both galaxies (top) and halo atmospheres (bottom), decom- 
posed into hot and cold contributions based on a compari- 
son of T max to T v ir iacc . Here, in contrast to earlier figures, we 
extend the "time window" over which accretion is counted 
to include a longer evolutionary history of these systems. 
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Figure 6. The total amount of gas smoothly accreted onto either galaxies or halo atmospheres by z = 2 decomposed into hot and 
cold contributions. Following our earlier discussion we make this separation based on some ratio of T max to T v i r acc , as shown. Unlike in 
previous figures, we include gas accreted over an "integrated" time window encompassing the evolution of these systems back to z = 6 
or as far as they have well-defined parents. 



We extend back in time as far as each halo is successfully 
tracked in our parent tree scheme, up to z = 6. Consequently 
Figure [6] shows the integrated contribution of smoothly ac- 
creted baryons to the z = 2 mass budgets of these systems. 
When compared to the accretion rates calculated only over 
the past 1 Gyr we find that the previous conclusions regard- 
ing the comparisons between the two codes and the balance 
between the hot and cold accretion modes remain qualita- 
tively the same. Since the two simulations have the same 
initial conditions, we naturally expect their properties to 
tend towards better agreement at sufficiently high redshift. 
In particular, comparing the total accreted mass onto galax- 
ies at halo masses of ~10 12 M Q we find that AREPO has a 
factor of ~10 higher mass accreted with T max > T v i r , and 
a factor of ~3 lower mass accreted with T max < T V - 1T when 
compared to GADGET. This represents a significant discrep- 
ancy in the gas accretion history integrated down to z — 2, 
which is directly related to the morphological and kinematic 
differences of the respective galaxy populations. 

A long standing difficulty of SPH has been its abil- 
ity to produce central objects with angular momentum suf- 
ficient to form rotationally supported disks with realistic 
scale lengths jVogclsber ger et al.| [2012). The relative differ- 
ences between the angular momentum content of cold gas 
accretion in GADGET when compared to hot gas accretion 
in AREPO explains why galaxy properties differ between 
these two approaches. AREPO galaxies are more disk-like 
and more rotationally supported, with systematically larger 
disk scale lengths ( |Torrey et al.|[20lT| ) in part because they 
are significantly built from the diffuse, hot halo gas, which 



has a large reservoir of angular momentum (see also Sijacki 
et al. 2012). In contrast, GADGET galaxies are more com- 
pact and less rotationally supported in part because they 
are being built mainly from the cold gas, which is clumpy 
due to numerical artifacts and which tends to artificially lose 
its angular momentum. That realistic disk morphologies are 
related to the angular momentum reservoir from which they 
form has also been found regardless of numerical technique, 



in the context of galactic fountain material (Brook et al 
2012 1 and cooling from hot coronae (Sales et al.|2012[ ). Here 



the differences in the hot halo itself may also play a role. 
As we discuss below, hot halo gas in AREPO has a smaller 
spatial extent and reaches somewhat higher central densities 
than in GADGET ( |Keres et al.|2012[ ), which leads to shorter 
cooling times and thus larger hot accretion rates onto cen- 
trally forming galaxies. 



3.4 Geometry of Accretion 

We find that the filamentary nature (geometry) of accreting 
gas at and around the virial radius of massive haloes is in 
good qualitative agreement between GADGET and AREPO. 
At high redshift (z — 2) filaments are a common feature in 
a majority of massive (M ^ lO 11,5 A/0) haloes. Prior to this 
section, we have included only smoothly accreted gas in our 
analysis and discussion. When visualizing the geometry of 
accretion, we show all gas elements, and do not restrict to 
only the smooth accretion mode as in previous plots. Fig- 
ure [7] shows a gas projection of a particular halo matched 
between the two simulations, where in the top two rows tem- 
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Figure 7. Density and temperature distribution of individual gas elements for a ~ 10 11,5 Mq halo at z = 2 in GADGET (left panels) 
and AREPO (right panels). The top panels show gas of all temperatures, while the middle panels show only gas with instantaneous 
temperature T < 10 5 K. Ticks are color coded by instantaneous temperature and with directions representative of the local velocity field. 
The bottom panels show the mass- weighted temperature projection for lines of sight through a larger cube of side length 5r v j r . The same 
dark matter halo in AREPO has a less extended hot component than its GADGET counterpart. In both cases locally overdense and 
overcool filamentary gas structures penetrate the hot halo at r v - lr and deliver gas to smaller radii. In GADGET these filaments become 
extremely collimated and cold at ~0.5r v ; r whereas in AREPO they tend to become more diffuse, dissipate, or experience significant 
heating. This last case can be seen at the top of the left filament which shows a thin, ~ 10 K streamer. 
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Figure 8. Density of individual gas elements for a ~ 10 10 - 7 A/ Q halo at z = 2 (T vir ~ 10 5 4 K) in GADGET (top panels) and AREPO 
(bottom panels) color coded by instantaneous temperature and with directions representative of the local velocity field. The left panels 
show gas of all temperatures, while the right panels show only gas with instantaneous temperature T < 10 5 K. This halo, as with a 
majority of lower mass systems, does not show the same prominence of filamentary gas inflows. 



perature is mapped to particle color and the local velocity 
field by tick direction. The top panels shows all temperature 
gas, while the middle panels show only gas with instanta- 
neous temperature below 10 5 K. In GADGET (left column) 
we observe that filamentary gas structures often coherently 
flow to r ^ 0.25r v i r and become progressively smaller in 
cross section. They can connect directly to central gas struc- 
tures at r ^ 0.1r v i r , although the contribution of these thin 
streams to the total accretion rate of gas onto the galaxy is 
unclear. Their temperature profiles remain well below that 
of the virialised halo gas and show little evidence for shock 
heating. In AREPO (right column) these same gas filaments 
are more diffuse and experience significant heating at com- 
parable radii. For instance, the two filaments from the top 
and left which seem to "disappear" in the right panel in 
fact heat to the temperature of the halo gas while remain- 
ing coherent. We discuss this important contribution to the 
accretion rates of hot gas below. 

Another significant difference between the two results 



is the distribution of the hot halo gas itself. The bottom 
row of Figure [7] shows the mass- weighted temperature pro- 
jection for the same halo in both AREPO and GADGET. In 
order to make a fair visual comparison we integrate each 
sightline by spatially distributing the gas mass of both SPH 
particles and Voronoi cells using the same SPH kernel ap- 
proach. We observe that hot halo gas in GADGET extends 
to larger radii than in AREPO. To quantify this effect, we 
measure the stacked radial temperature profile of all haloes 
in the box with lO n ' 5 M < M ha lo < 10 12 ' M Q . At the virial 
radius, the gas temperature in AREPO haloes is a factor of 
~ 2/3 lower than in GADGET - the mean gas temperature 
at r v ir in GADGET is obtained at ~ 0.7r vir in AREPO. As 
we discuss in the context of accretion rates, the smaller hot 
gas reservoir in AREPO haloes results from more efficient 
cooling out of the halo, while in the case of GADGET this 
cooling is suppressed due to spurious numerical heating. 

An independent point which has led to confusion is the 
nature of accretion in lower mass haloes. Figure [8] shows 
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Figure 9. An all-sky Mollweide projection of temperature, overdensity, and radial mass flux on the surface of the half virial sphere 
centered on one halo with log(M) ~ 11.8. All gravitationally bound substructures have been removed. The clustered regions of negative 
radial mass flux represent cross sectional views of gas filaments. They are associated with mass overdensities and cold temperatures with 
respect to the shell averages. The larger arcing features in the map are infalling sheets of gas with large angular extent on the sky. 



an example of such a system with Mhaio — 10 1 Mq at 
redshift two. In this case there are no prominent filaments, 
nor do we expect any ( Katz et al.||2003 Keres et al.|2009 |. 
The virial temperature of this system is ~10° * K and would 
lie at the low mass end of the transition mass region from 
cold to hot mode dominated. Given the usual definition of 
logT c = 5.5, gas which had accreted and shock heated to the 



virial temperature would be considered cold mode accretion. 
In contrast, under a definition based on some fraction of T v i r 
a large fraction of that same gas would be considered hot 
mode. Indeed, for this mass range in Figure[5]we expect that 
approximately half of the halo gas has T max > T v i r and half 
has T max < Tvi r . It is then critical to differentiate between 
gas which shocks to T v i r ^ T c and gas which accretes "cold" 
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Figure 10. Comparison of strictly cold gas with instantaneous temperature below 10 5 K for eight matched systems at redshift two. 
The top eight panels show systems from the GADGET simulation while the bottom eight panels show the identical systems in AREPO. 
We commonly find examples of predominantly spherical accretion in low halo masses and differences in filamentary gas structures for 
larger halo masses, as discussed in the text. 



into low mass haloes with short cooling times. This is par- 
ticularly the case since the signature of both mechanisms is 
similar in maximum past temperature, making any separa- 
tion based on temperature alone difficult if not impossible. 

In order to better compare the properties of filaments 
we construct maps of various fluid quantities which re- 
veal the cross-sectional profiles of these structures. Figure 
[9] shows gas temperature, density and radial mass flux cal- 
culated on the surface of a sphere centered on one particular 



halo, the log(M) = 11.8 example from Figure 10 and with 
a radius equal to half the virial radius. We use a simple 
mass-weighted tophat kernel to interpolate these quantities 



onto equal area pixels using the HEALP1X scheme (Gorski 
ct al. 2005), and display the result using a Mollweide all-sky 
projection. 

Pronounced differences appear in images of projected 
gas density. We find that the radial mass flux inward at 
these radii is dominated in the SPH halos by a large num- 
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ber of thin filaments and gas blobs with progressively smaller 
cross section. The moving mesh case shows no evidence of 
gas structures of comparably small angular extent. The den- 
sity profiles of the filaments are also not as centrally concen- 
trated and have larger angular extent. Although filaments 
reach roughly the same minimum temperature, the transi- 
tion to the hot 10 5 ' 5 — 10 6 K gas in the halo is more gradual 
in AREPO. That is, filaments have less well-defined bound- 
aries, indicative of increased mixing, in the moving mesh 
calculation, though they appear to be roughly at the same 
level of pressure equilibrium. Although the angular covering 
factor of infalling material is slightly larger in AREPO, it is 
interesting that for this particular halo, where we have ac- 
cretion onto the central galaxy dominated by T max > T v i riacc 
gas, the mass flow is still highly aspherical. That these fil- 
aments of hot gas, originating from large scale features in 
the IGM, contribute to the hot gas accretion rate distinctly 
from classical cooling flows out of hot halo gas is the third 
key point of this work. The presence of this coherent hot ac- 
cretion implies that a simple cooling rate calculation based 
on halo properties cannot accurately describe high redshift 
accretion in massive haloes. This may have important con- 
sequence for semi-analytical models (SAMs) which assume 
radiative cooling from a quasi-static atmosphere. 

Comparing large ~ f Mpc scale gas and dark matter 
density fields, the cold ~ 10 4 K filaments which are rela- 
tively narrow within the virial radius can clearly be iden- 
tified as extensions of filaments of larger diameter outside 
the halo. Given that GADGET and AREPO share the same 
treatment of dark matter and the same gravity solver, and 



merging substructures, both resolved and unresolved (Keres 



show good agreement in IGM properties outside haloes ( Bird 



|et al. ||2012p , it is unsurprising that at the virial radius the 
differences we observe above minimise and the two codes 
show similar gas properties. 

We show eight further examples of strictly cold gas 
below 10 5 K in matched systems spanning the mass range 
10.7 < log(M ha i o [M ]) < 12.0 in Figure [To| illustrating that 
the similarities and differences we describe are commonly 
found in the simulated cosmological volume]^] Haloes in both 
simulations support approximately the same frequency of fil- 
amentary structures over any given range of halo mass. The 
exception appears to be filaments that are tidal in nature, 
which are much more pronounced - longer and more con- 
centrated - in AREPO. The cold gas distribution in the most 
massive GADGET haloes is dominated by the large blob pop- 
ulations, which are absent in the corresponding AREPO sys- 
tems. 



3.5 The Clumpy Contribution 

We seek to understand the accretion of material onto 
haloes and galaxies originating directly from the intergalac- 
tic medium. If we include all material, the resulting cold 
fractions and cold mode accretion rates are only upper lim- 
its, due to some nonzero contribution from minor mergers. 
Past numerical simulations run with SPH and focusing on 
cold mode accretion have specifically addressed this issue of 



2 A catalog of similar visual comparisons for all objects with 
M > 10 10 - 75 M Q at z = {0,1,2,3} is available online at 
www. cf a. harvard . edu/itc/research/movingmeshcosmology/. 



et al.|[2005] 2009| |Brooks et a"L] [20091 |van de Voort et al 



2011 1. The exact approach in separating out the merger con- 

and in this study we 



tribution differs. InjBrooks et al 



(2009) 



remove material that belonged to any galaxy halo other than 
the main galaxy under consideration, while in Kcrcs et al. 
( 2009 ) for instance only the central galaxy and not its asso- 
ciated circumgalactic (or halo) material is removed. Other 
studies have neglected this distinction and included the ac- 



cretion of dwarf satellites as part of a cold mode ( Shen et al 



2012|. Although simulations using the AMR technique note 



the presence of substructure embedded within filaments, 
they also neglect to make this separation ( Dekel et al.|2009 



Agertz et al. 



2009 1 though due to technical constraints. One 



option as in |Ocvirk et al. ( 2008 1 is to use an instantaneous 
criterion such as gas density to identify a clumpy compo- 
nent. However, without some additional tool for tracing gas 
properties back in time (tracer particles) the contribution 
of substructures and associated material cannot be unam- 
biguously accounted for, and a measurement of the cold gas 
accretion from the IGM represents an upper limit. 

Applying the three definitions for "smooth" , "clumpy" , 
and "stripped" from Section |2.5| we decompose the accre- 
tion rates separately for hot and cold accretion in Figure [TT) 
where we here use 1.0 x T v i r ,acc as the separating criterion. 
As previously discussed, our definition of clumpy includes all 
material bound to satellites - both galaxies and their associ- 
ated circumgalactic medium. We focus first on cold material 
accreted onto central galaxies, for which the clumpy mode 
is the ~ 60% dominant contributor at all halo masses in 
both GADGET and AREPO. Smooth accretion accounts for 
only ~ 10% of the cold mass accumulation at the massive 
end, where most of the cold accretion onto central galaxies 
arises from merging substructure. For the central galaxies 
of the lowest mass systems we consider with ~ 10 10 Mq the 
smooth contribution is slightly larger (~ 20%). This frac- 
tion is relatively converged with resolution; however, any 
estimate of the merger contribution for the lowest mass sys- 
tems may be underestimated due to numerical resolution. 
For all halo masses the stripped component is comparable 
to the smooth component. The accretion of hot material 
onto central galaxies shows similar behavior, except that 
the contribution from clumpy mode gas scales up with in- 
creasing halo mass, while the contribution of smooth gas 
scales down. The only notable difference between the two 
codes is for material accreted hot onto galaxies hosted in 
massive haloes with Mhaio 10 11 A/ Q for which AREPO 
shows a lower clumpy fraction (40% vs 70%) and a higher 
smooth fraction (40% vs 20%), consistent with the picture 
that smooth gas is more efficiently heated in AREPO. 

The importance in separating different gas origins in 
the simulations can be seen by reference to observations of 
our local environment, in which the Magellanic Stream is a 
large reservoir of cold, neutral gas thought to arise from the 
tidal interaction of dwarf satellites ( jBesla et al.|2010{|2012[ ). 
Its significant spatial extent through the halo of the Milky 
Way does not arise from a gas filament of direct cosmologi- 
cal origin, although it might easily be characterised as such. 
Similar features in external galaxies will be even more dif- 
ficult to characterise. Although separating the contribution 
of minor mergers and accreted tidal debris from gas infalling 
directly from the IGM presents a technical challenge for sim- 
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Figure 11. The total accretion rate separated into contributions from our three different modes: smooth (dotted), clumpy (dashed), 
and stripped (dot-dashed). 



ulations, it is necessary to truly identify the origin of the gas 
that ends up forming galaxies, without mixing up contribu- 
tions from disparate physical origins. 



3.6 Resolution Convergence 

In this last section we address some of the numerical is- 
sues related to our moving mesh calculation. In Figure [12] 
we show the convergence of the gas accretion rate onto cen- 
tral galaxies for our three resolution levels. The total (the 
sum of the top and bottom panels for any given column) 
is well converged across the given halo mass range, and all 
resolution trends are qualitatively similar regardless of the 
separating line drawn between hot and cold. In general, |Vo-| 



gelsberger et al. (20121 showed the increase in high redshift 



global star formation rates due to better resolving dense gas 
in low mass haloes. We expect that smooth accretion rates 
will then decrease as resolution increases, with a correspond- 
ing increase in either clumpy or stripped accretion, as more 
gas mass collapses into substructures at early times. Galax- 
ies in AREPO exhibit this behavior for material accreted 
cold, i.e. directly from the IGM. In contrast, the accretion 
rate of hot material from the halo increases strongly with 
resolution for halo masses which are poorly resolved in our 
simulations. This enhanced cooling from the halo occurs due 
to a combination of the changing temperature distribution 
of hot halo gas, together with inadequately resolved mixing 
in two particular regimes, as we discuss below. 



First, Keres et al. (20121 shows that the radial temper- 



ature profiles of AREPO haloes become steeper and reach 
higher maximum temperatures with increasing resolution. 
We verify this and furthermore find that it is particularly 



problematic for poorly resolved systems where the number 
of gas cells in the entire halo is of order a few thousand or 
fewer. For more well-resolved haloes (Mhaio 10 Mq) the 
temperature profiles are converged between our medium and 
high resolution runs. This occurs in part because the halo 
material more fully populates the high temperature tail of 
the thermal distribution of the virialised gas. 

The strong mixing between the rotating galactic disk 
and the halo gas identified by |Sijacki et al.| ( |2012| ) is also 
artificially enhanced at low resolution where this interface is 
poorly resolved. This shortens the cooling times of hot halo 
gas near the disk and drags down the temperature profile, 
which ultimately results in accreting material heating to a 
lower peak temperature. The potential for overmixing at low 
resolution is a deficiency inherent to grid codes, including 
AREPO. In this case the moving mesh nature of the scheme 
acts to minimise the effect, though it cannot be suppressed 
entirely as in SPH codes. The same issue arising at the disk- 
halo interface also occurs in a second regime, with infalling 
substructures which penetrate into the halo and lower the 
cooling times of the hot gas due to overmixing. This is in 
contrast to SPH, where a relatively poor treatment of fluid 
instabilities prevents efficient stripping of gas from infalling 
satellites (Sija cki et al.|2012[ ). 

Galactic disks in the low mass haloes are also poorly 



resolved in our lowest resolution simulations (Keres et al 



2012 1, which may extend them beyond our 0.15r v j r cut, par- 



ticularly in AREPO where they are larger and more extended 
to start with. The result is an underestimation of accretion 
rates, while at higher resolution these disks are more com- 
pact and the accretion is largely recovered. This effect may 
preferentially influence gas with larger T max which is found 
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Figure 12. Convergence of the accretion rate onto galaxies with resolution. The vertical black lines show the halo masses at which 
^b-Mhalo/ m targct - 5 X 10 4 for the 512 3 and 256 3 runs, a rough requirement on the number of gas mass elements for which gas accretion 
physics begins to be adequately resolved and converged in AREPO. 



at systematically larger radii and lower density in galactic 
disks. 



4 DISCUSSION AND CONCLUSIONS 

We are motivated to define a cold accretion mode based on 
the markedly distinct channel of cold filamentary streams 
penetrating the hot atmospheres of massive haloes. Given 
the existence of two distinct modes, how can we distinguish 
between the two? We observe that this form of gas accretion 
enables gas to infall onto galaxies without shock heating to 
an appreciable fraction of the virial temperature. This mo- 
tivates a selection criterion based on virial shocking, or lack 
thereof. This is not a new idea. Indeed 



Keres et al. ( 2005 1 



explored this option but concluded that the distinction be- 
tween hot and cold was less clear than with a constant tem- 
perature cut. van de Voort et al. (20111 likewise compared 
hot fractions obtained using either a threshold logT c = 5.5 
or various fractions of T max /r v i r . That work concluded, as 
we do, that the hot fraction depends strongly on the choice of 



criterion. Furthermore, van de Voort et al. (20111 acknowl- 



edges that gas accreted in their "cold mode" may have in 
fact gone through a virial shock, indicating that it is critical 
to investigate the gas selection made with a global constant 
temperature threshold over the range of halo masses being 
considered. 

Furthermore, even comparing the temperature history 
of gas elements to the virial temperatures of their parent 
haloes is an inexact probe of whether or not gas shock heats. 
The picture that gas shock heats up to T v i r at r v i r is an over 



simplification of a virialisation process that may involve mul- 
tiple shocks at increasingly higher temperatures as infalling 
material encounters the temperature gradient of halo gas. 
Indeed, shock heating may not even be the most funda- 
mental question as to how galaxies obtain their gas. This 
might be better answered by considering the cooling time 
of gas ( |White fc Rees|p78l | White fc Frenk]|T99Tj ) , and so 
the time-scale for collapse onto a forming galaxy. The radial 
stalling associated with incorporation into the hot halo gas 
may be more useful than the gas temperature itself. In the 
same manner, the angular momentum history of infalling 
gas may provide a more physically motivated segregation of 
different accretion modes. We have found that the interpre- 
tation of the angular momentum, for instance, is not nearly 
as clear as with maximum past temperature and we suggest 
this is the primary reason why T max remains the predomi- 
nant technique used to study cold mode accretion. 

Based on additional temperature cuts, recent simula- 
tions have also identified a third "warm" mo de of gas accre- 
tion ( Joung et al.|2012" 



Murante et 



in the 10 5 K < 

T < 10° K regime and selected by constant temperature 



al.|2012 l 



bounds over either instantaneous temperature (in Joung 
et al. (2012 I using AMR) or maximum past temperature (in 
Murante et al.| ( |2012[ ) using SPH). In the case of |Murante] 
et al. (20121 this prominent mode is attributed to effective 



SN thermal feedback. Unfortunately, since these studies are 
focused on low redshift and include galactic outflows we can- 
not make a direct comparison to our results. We have not 
therefore investigated this particular temperature selection 
criterion, although a quantitative comparison with AMR re- 
sults and an exploration of the impact of feedback are cru- 
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cial and natural directions for future work. The commonly 
perceived agreement between SPH and AMR simulations, 
of cold mode accretion as a dominant contributor to galaxy 
formation, is based on a mixed analysis of instantaneous gas 
properties and qualitative comparisons. An as of yet unper- 
formed comparison of the thermal history of accreted gas 
(using tracer particles) is required to demonstrate quantita- 
tive agreement between these two methods. 



4.1 Numerical Uncertainties and Included Physics 

Particularly when studying gas accretion, a perceived 
strength of SPH, embodied in codes such as GADGET, is 
what should accurately be termed its pseudo-Lagrangian 



et al. ( 2009 1. However, this agreement is a coincidental result 



nature ( Vogelsberger et al. 2012 1. The consequence of this 



behavior is the apparent ease by which the thermodynamic 
properties of individual mass elements can be traced through 
time. However, consider the example of mass represented 
by a single SPH particle in a shearing flow, as illustrated 



in Fig. 20 of Vogelsberger et al. (20121. This mass is ini- 



tially drawn from a volume centered on the SPH particle 
with radius of order the local smoothing length. Shearing 
of this volume deforms its shape and the associated mass 
should follow. However, the nature of SPH is such that this 
mass is forced to remain bound to its SPH particle. The 
consequence is that the evolution of the fluid is not fully 
consistent with the equations of motion and cannot then be 
truly Lagrangian. In this regard AREPO is also not strictly 
Lagrangian since the gas cells are not allowed to become ar- 
bitrarily distorted. However, the fundamental difference is 
that mass exchange between cells occurs in a manner con- 
sistent with the equation of mass conservation and we can 
therefore term the scheme quasi-Lagrangian, because unlike 
in SPH the solution is faithful to the underlying equations 
of motion. 

It is also important to note the differences arising due 



solely to different formulations of SPH. Keres et al. ( 2005 1 
used TreeSPH ( Hernquist fc Katz|1989 1 and found a strongly 
bimodal histogram of the maximum past temperature of ac- 
creted gas, and comparable contributions of the hot and 



cold modes. Keres et al. (20091 used GADGET-2 and con- 



cluded that at z ^ 2 the accretion is dominated by the cold 
accretion mode at all halo masses, with little evidence for 
bimodality in the same T max histogram. The original moti- 
vation for a constant threshold temperature used to separate 
hot and cold modes is in fact the location of the temperature 
minimum resulting from this bimodality. Yet, this same con- 
stant value has been used in virtually all simulation based 
investigations of cold mode accretion since. 



Yoshida et al. (2002) compared these two SPH formu- 



lations and concluded that the TreeSPH technique leads to a 
significant over-cooling stemming from its geometrical sym- 
metrization of hydrodynamical forces. In the context of cool- 
ing from halo gas this problem manifests as increased hot 



mode accretion. Accordingly, Keres et al. ( 2009 1 attributes 



the much higher hot mode accretion rates in their earlier 
simulations to numerical deficiencies in the older SPH for- 
mulation. The larger relative contribution of high T max gas 
and the dominant contribution of hot gas accretion in the 
most massive haloes found in |Keres et al.|p005[ ) is in fact in 
better agreement with this work than that found in |Keres] 



of the aforementioned numerical inaccuracies in TreeSPH. 

It is not altogether surprising, then, that our moving 
mesh results with AREPO are in tension with previous stud- 
ies, particularly SPH simulations made using the standard 
formulation found in GADGET. Significant uncertainties re- 
main in the numerical simulation of cosmological volumes 
where the inclusion of a diverse number of physical processes 
over a large dynamic range in spatial and temporal scales 
is both computationally intractable and yet phenomenologi- 
cally required. What is surprising is the extent to which dif- 
ferences arise solely due to the numerical scheme employed 
to solve the hydrodynamics, and not to differences in sub- 
grid prescriptions, star formation recipes, feedback imple- 
mentations, and the like. The simulations presented in this 
work implement a basic set of physical processes in addi- 
tion to gravity and hydrodynamics. We leverage this rela- 
tive simplicity to make a clean comparison between the SPH 
and moving mesh techniques. However, it is probable that 
strong feedback from stars or active galactic nuclei could sig- 
nificantly alter the nature of gas inflow from the IGM, and 
this remains an open question for future work. The mov- 
ing mesh implementations of ideal magnetohydrodynamics 



(Pakmor et al. 20111 and physical viscosity (Munoz et al. 
2012[ ) in AREPO will likewise enable us to investigate these 
potentially important effects for cosmological gas accretion. 



4.2 Conclusions 

We conclude with two main points. Firstly, that the constant 
temperature threshold commonly used to study cosmolog- 
ical gas accretion is only reasonable above some minimum 
halo mass, and application in an overly broad context bi- 
ases conclusions regarding the relative importance of hot 
or cold gas accretion. Secondly, that numerical deficiencies 
inherent to smoothed particle hydrodynamics (SPH) simu- 
lations non-trivially modify the relative contributions of hot 
and cold mode accretion, under any definition. An identical 
analysis of gas accretion utilizing our Monte Carlo tracer 
scheme with the new moving mesh code AREPO demon- 
strates significant physical differences in the thermal history 
of accreted material. We summarise our primary results as: 

• In agreement with previous work, GADGET simulations 
imply that at z — 2 only a small fraction of gas in centrally 
forming galaxies of massive haloes above 10 11 Mq heats to 
an appreciable fraction of the virial temperature during ac- 
cretion. In AREPO we find a decrease in the accretion rate of 
cold gas, by a factor of ~ 2 at Mhaio — 10 11 Mq, as well as a 
significantly enhanced accretion rate of hot gas, by an order 
of magnitude at the same halo mass. These discrepancies in- 
crease for more massive systems. We attribute the decrease 
in the cold accretion rate primarily to the large population 
of numerical "blobs" which efficiently deliver cold gas to 
central galaxies in our SPH simulations, but are completely 
absent in our moving mesh calculation. The increase in the 
hot accretion rate is dominated by significantly more effi- 
cient cooling from halo gas in AREPO, where spurious heat- 
ing from the dissipation of turbulent energy on large scales 
prevents the same behavior in GADGET. 

• We argue that comparison of the maximum past tem- 
perature T max of a gas element to a fixed temperature 
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threshold T c makes physical sense only for haloes with 
Tvir >• T c . For lower mass systems the past temperature 
history should instead be compared to some fraction of the 
virial temperature of the dark matter halo. However, at suf- 
ficiently low halo masses, when the virial temperature be- 
comes comparable to the gas temperature in the IGM, nei- 
ther a constant temperature threshold nor one scaled with 
the virial temperature are sufficient to probe gas shock heat- 
ing and virialisation. 

• We observe that the "transition mass" from cold to 
hot dominated accretion which has been reported between 
1O 11_11 ' 5 M0 is a consequence of the constant temperature 
criterion. When comparing the thermal history of gas in- 
stead to a fraction of T v i r we find no sharp transition, only 
a gradual decline from / co id — 20% to ~ 0% over the mass 
range from M ha i G = 1O 1O M to 1O 12 M . 

• The filamentary geometry of accreting gas near the 
virial radius is a common feature of massive haloes above 
M ~ 10 11 ' 5 at high redshift (z — 2). Although characterised 
by the same large scale morphology, filamentary gas struc- 
tures in GADGET tend to either remain cold and flow co- 
herently to small radii within a halo, or artificially fragment 
and form a large number of "blobs" which are purely nu- 
merical in origin. In contrast, filaments in AREPO simula- 
tions are more diffuse and experience significant heating at 
comparable radii. The geometry and angular covering fac- 
tor of material accreted with T max > T v ir,acc indicates that 
coherent, filamentary flows associated with large scale IGM 
filaments contribute significantly to hot accretion rates in 
these massive systems. 
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